Spin relaxation in diluted magnetic semiconductor quantum dots 
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Electron spin relaxation induced by phonon-mediated s-d exchange interaction in a II- VI diluted 
magnetic semiconductor quantum dot is investigated theoretically. The electron-acoustic phonon 
interaction due to piezoelectric coupling and deformation potential is included. The resulting spin 
lifetime is typically on the order of microseconds. The effectiveness of the phonon-mediated spin- 
flip mechanism increases with increasing Mn concentration, electron spin splitting, vertical confining 
strength and lateral diameter, while it shows non-monotonic dependence on the magnetic field and 
temperature. An interesting finding is that the spin relaxation in a small quantum dot is suppressed 
for strong magnetic field and low Mn concentration at low temperature. 
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I. INTRODUCTION 

The spin of the electron in low-dimensional semicon- 
ductor structures has received intense interest in re- 
cent years due to its potential applications in spin- 
tronic devices and quantum information processing 
technologies^*^ To improve the performance of such 
devices, the decoherence of the electron spin due to 
coupling to environmental degrees of freedom should 
be minimized. Theoretica]i^*S*S*i and experimental 
investigations^i^ have shown that the electron spin could 
have a extremely long relaxation time in nonmagnetic 
semiconductor quantum dots (QD's), compared with 
that in the bulk or quantum wells. Theoretical works pro- 
posed that the diluted magnetic semiconductor (DMS) 
QD's can be used as spin aligners, spin memories as well 
as spin qubits,i° DMS QD's offer us a new flexibility 
in manipulating carrier spins, since the spin properties 
can be strongly influenced by applying an external mag- 
netic field or varying the temperature.^^ Various relax- 
ation mechanisms of the electron spin come from differ- 
ent coupling to the environment, i.e., magnetic impurity, 
nuclear spin, and spin-orbit interaction. It is important 
to identify the dominant mechanism of spin relaxation 
for a particular system. Previous theoretica U-f - " and 
experimentaliiii^ works on spin relaxation in DMS quan- 
tum wells have indicated that the s-d exchange interac- 
tion between band electrons and localized spins of mag- 
netic ions is the dominant spin-flip mechanism, leading to 
electron spin lifetime of the order of picoseconds. How- 
ever, to the best of our knowledge, there is no study on 
the spin relaxation induced by the s-d exchange interac- 
tion in DMS QD's. Detailed theoretical and experimental 
investigations are necessary to gain physical insight into 
the spin relaxation process in such structures. 

In this paper, we investigate theoretically the spin re- 
laxation of the lowest Zeeman doublet in vertical II- VI 
DMS QD's. The spin-flip scattering caused by the acous- 
tic phonon-mediated s-d exchange interaction between 
the conduction electron and Mn ions is considered. The 
electron-acoustic phonon coupling includes the piezoelec- 
tric and deformation potential interactions. Since the 



first-order spin-flip process through direct scattering by 
the Mn ions is generally blocked by the energy-matching 
condition, we consider the second-order process involv- 
ing the emission or absorption of a phonon. Our calcula- 
tion shows that this phonon-mediated spin-flip scattering 
leads to electron spin lifetime typically of the order of mi- 
croseconds. The effectiveness of this mechanism increases 
signiflcantly with increasing Mn concentration, electron 
spin splitting, vertical confining strength and lateral di- 
ameter, while it shows non-monotonic dependence on the 
magnetic field and temperature, due to the competing ef- 
fect of the electron spin splitting, the phonon number and 
the correlation function of the Mn ions. It is interesting 
to notice that the spin relaxation of electrons in the low- 
est Zeeman doublet is suppressed in a strong magnetic 
field at low temperature for II- VI DMS QD's with low 
Mn concentration. 

The rest of this paper is organized as follows: the theo- 
retical model and formula of the spin-flip scattering rate 
(SFR) induced by phonon-mediated s-d exchange inter- 
action are derived in sec. II. Numerical results and dis- 
cussions for the SFR as a function of magnetic field, as 
well as its dependence on the Mn concentration, QD size 
and temperature are given in sec. Ill and we give a brief 
conclusion in sec. IV. 



II. THEORY 

We consider II-VI DMS QD's subjected to a per- 
pendicular magnetic field. Assuming an infinite deep 
well along the growth direction (the z axis) and a in- 
plane parabolic confining potential, the electron wave 
function can be written as ^j{v) = x{z)4'{Pt where 
xiz) = y^2/ zq sin(7rz/zo) is the ground state wave func- 
tion along the z axis (we have assumed that the vertical 
confinement is strong and only the lowest energy level is 
relevant), zq is the width of the well, and (j}(p, cp) is the in- 
plane wave function determined by the two-dimensional 
Hamiltonian H = Hq + Hg-d + He-p- The first term 

H, = ^^^^ + \^*^Ip' + \9*P^BBa^ (1) 
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is the electron Hamiltonian in the external magnetic field 
and parabolic potential. Here m* is the electron effective 
mass, A = {—By/2,Bx/2,0) is the vector potential, llIq 
characterizes the lateral confinement strength, g* is the 
intrinsic electron g- factor, and az is the z-component of 
the Pauli matrices. The second term 



J(r - R,)s • 



(2) 



Now the total Hamiltonian is divided into two parts, 
= + where 



— ^ + -TO CJ p + -i^cl^z + T^CT^A^, (5) 
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describes the s-d exchange interaction between the elec- 
tron (s) and the localized Mn ion (Si), where J(r) is the 
s-d coupling integral, and the summation runs over all 
the Mn sites. The last tcrm^ 
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describes the interaction between the electron and acous- 
tic phonon, where h^vih^^) is the annihilation (creation) 
operator of the bulk phonon mode with wave vector q 
and branch v. 

The s-d exchange term is divided into a mean-field 
part and a fluctuating part, Hs-d = H^_^ + Vs-dj where 

H°_^ = azAsd/2, Vs-d = - E.Jir - R0(-s^+^5^^ + 
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Sf±S^, Asd 



—Noax (Sz) is the exchange splitting, Nq is the number 
of unit cells per unit volume, a — {cf)c | J(r)| 0c) 
is the unit cell volume, (pc is the Bloch function at the 
bottom of the conduction band) is the s-d exchange cou- 
pling constant, x is the fractional occupation factor of 
the cation sites by the Mn ions. 
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is the thermal average of the Mn spin, with S ~ 5/2 the 
Mn 3d^ spin, Bs{x) the Brillouin function, and 5*0, Tq 
phenomenological parameters accounting for the antifer- 
romagnetic superexchange between neighboring Mn ions. 



V = Vs-d + He-p. Here u = uo^ ojI/A, Uc = eB/m* 
is the cyclotron frequency, Lz is the z-component of the 
orbital angular momentum, = g* -\- Ag^ is the 
total Zeeman splitting of the electron. 

In order to obtain the SFR induced by the s-d ex- 
change and electron-acoustic phonon interaction, we con- 
sider the whole system (including the electron, Mn ions 
and the phonon bath) transits from an initial state \i) = 
\la; M ; N) to aU possible final states |/) = M'; N') in 
which the electron spin is reversed. Here \la) is the elec- 
tron eigenstate {I stands for the orbital quantum number 
(n, m), see the Appendix for details), a — ± (ct = =f) de- 
note spin-up (spin-down) and spin-down (spin-up) state, 
respectively. \M) = jMi^, M2z, • • • ) is the eigenstate of 
the Mn ions and |A^) — Y\^^ W^y) denotes the phonon 
state. The transition rate is averaged over the random 
positions and the initial states of the Mn ions, as well as 
the initial states of the phonon system to give the SFR 
of the electron from the initial state \1(t) to the final state 
|/'(7), denoted as Wi'gja- 

Since the spin-flip process of electron is always accom- 
panied by the flip of a Mn spin due to the conservation of 
the total angular momentum (see Eq. l|2Il), we introduce 
the renormalizcd electron energy Ei± = ejit Ao/2, where 
El is the orbital eigenenergy of Hq (see the Appendix for 
details) and Aq = Agd — A^ is the (renormalized) elec- 
tron spin splitting, with A^ — {gMn — g*)tJ-BB. Based 
on second-order perturbation theory, the transition am- 
plitude between |i) and |/) is given by 
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El'g- — Elj^a 

{I'a; M'; N'\He-p\lia; M'; N) {ha; M'; N\Vs^d\l(J; M; N) 
Eia - Evs 

I 



(6) 



In the first term, the electron first hops from the initial 
state I^ct) to an virtual state with the same spin \h<^) 
through the interaction with a phonon, then it makes a 
spin-fiip transition to the final state \l'a) through the s-d 
exchange interaction with one Mn ion. The second term 
describes the process that the electron is first scattered 
to an opposite-spin virtual state \h&) through the s-d 



exchange interaction with one Mn ion, then it transits to 
the final state \l'a) via the assistance of a phonon. 



The scattering rate of the whole system is obtained 
from the Fermi golden rule Wfi = (211 jK) |r/i|^(5(£'/ — 
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Ei), and the SFR for the electron system is given by 



,i+ = -xiNoafG- 
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1 + sign(A/, 



1 - sign(Ai, 



Kii' 
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where 0'+ = {S'--^S'-+^) , G+- = {S'-+^S'-~^) are cor- 
relation functions of the Mn ions, S-f^) - 5^ ± iSy, 
^iv — El — Ell + Ao is the electron energy detuning, 
n{E) = [exp(— £^/(fcsr)) — is the phonon distribu- 
tion function, sign(a;) = 1 for x > and —1 for a; < 0. 
The kernel Kui is given by 



h-.h 
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where 



is the dimensionless overlap integral, and 



n J d^-R (Zi|R) (/2IR) (R|/3> (R\U) (10) 
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is related to the spin-conserved phonon-induced tran- 
sition rate. The explicit expressions for Si^^i^i^i^ and 
riiZ2Z3i4 are given in the Appendix. 

The spin lifetime Tia- of a given energy level \la) is given 

by 



Ti - ^— ^ 



(12) 



i.e., the sum of the spin-flip scattering rates from \la) to 
all opposite-spin final states ll'a). 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

From Eq. @, we notice that the kernel Kui and the 
SFR I/tjct (see Eq. (221) diverges when the energy of the 
intermediate state coincides with the initial or final state. 
To remove this divergence, we take into account the fi- 
nite lifetime of the intermediate level and add a small 
broadening parameter (an order-of-magnitude estimate is 
0.1me\4iS*i2'i£) to the energy of the intermediate state to 
convert this divergence into a resonance near the degen- 
eracy point lis This broadening parameter is not crucial 
for our calculation and would not change the qualitative 
behavior of the SFR. 

We consider Cdi_a;Mn2;Te QD's and use the following 
parameters in our numerical calculations: m*=0.096mo 



(mo is the free electron mass), 5*=— 1.6, CdTe mass 
density p=5.86 g/cm'^, lattice constant a=0.6481 nm, 
5Af„=2, S=5/2, Noa^220 meV, /ii4=0.394x 10^ V/m, 
sound velocity C/=3083 m/s, C(=1847 m/s. The lateral 
confining strength of the QD is characterized by the lat- 
eral diameter d = 2y^h/ {m*uj(j). The dependence of ^o. 
To on the Mn concentration x is obtained from Ref. 20. 

Considering the electron occupies the lowest spin-up 
and spin-down levels (n ~ 0,m — 0, ±), i.e., the lowest 
Zeeman doublet in the DMS QD, the SFR's can be cal- 
culated for the spin-up and spin-down states, which are 
denoted by 1/t+ and 1 /r_ , respectively. 

Since the spin-flip transitions to excited orbital levels 
are energetically unfavorable, the SFR of the lowest Zee- 
man doublet is dominated by the transition between the 
doublet, i.e., Woo+,00- and Wqo-,qo+, such that the elec- 
tron spin splitting Aq — Aoo,oo and the kernel i^Too.oo 
are important quantities for 1/t± (see Eq. (UJ and Eq. 
0). From Eq. (jnj, we note the contribution to Kqq^qq 
comes mainly from the intermediate level whose energy 
is the closest to the initial (or final) state (i.e., the orbital 
state (0, —1) in most cases, see the Appendix for details). 
That is, the term containing S'oo, 0-1,00, o-iFoo.o-i, 00,0-1 
and ^o-i, 0-1,00, ooFoo, 0-1,0-1, 00 (from the Appendix, we 
see they are equal to each other, so both terms are 
denoted as ST for short) in Eq. is the dominant 
contribution to -ftToo.oo- Additional contributions to the 
SFR 1/t± are the correlation function G '", G"' , the 
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FIG. 1: (a) Magnetic field dependence of tlie Mn correlator 
G ^ (solid curve) and G"* (dashed curve) at T=l K. (b) 
Phonon absorption (solid line) and emission (dashed line) fac- 
tor, (c) and (d): Piezoelectric coupling (solid curves) and de- 
formation potential (dashed curves) contribution to the prod- 
uct ST as a function of zq and |Ao|, respectively. 



phonon emission factor n(|Ao|) -I- 1 or absorption factor 
n(|Ao|). Therefore, in Fig. 1(a), (b), (c), (d), we plot 
schematically the correlation function G and , 
the phonon emission (absorption) factor and the prod- 
uct ST as a function of magnetic field, phonon energy, 
zq and |Ao| , respectively. In Fig. 1(a), we see that 
= S{S + 1) - {Sf) - \{S^)\ decreases monotoni- 
cally to zero while G-+ = S{S + 1)- {Sl) + \ (S,) \ shows 
a peak and approaches a constant value with increasing 
magnetic field. Physically, this is because G^ (G ^) is 
related to the transition of an electron from a spin-down 
(spin-up) initial state to a spin-up (spin-down) final state 
(see Eq. Q and Eq. ((HJ). Due to the conservation of 
the total angular momentum, the ^-component of the 
Mn spin (Sz) should decrease (increase) by one in this 
process. In a strong magnetic field, however, all the Mn 
spins are polarized antiparallel the magnetic field (i.e., 
{Sz) ~ —5/2), thus the correlation function G~^~ tends 
to vanish and the spin-flip process of the spin-down state 
is suppressed. The decrease of ST with increasing zq (ap- 
proximately ST cx 1/zo) E^nd the peak behaviors of ST 
as a function of | Aq | can be appreciated from Eq. (jA.2p , 
(jA.3|l and (jA.4|l in the Appendix. Note in the region 
where |Ao| is small, the dependence of ST on |Ao| is in 
agreement with Ref. 4. A peculiar feature is ST vanishes 
when the electron spin splitting Aq approaches zero, due 
to the vanishing energy and, as a result, the vanishing 
density of states of the involved phonon. 

Next, we shall investigate 1/t± as a function of mag- 
netic field for different temperatures, Mn concentrations 
and lateral diameters. The effect of vertical confining 
length Zq, lateral diameter d and temperature on 1/t± is 
also presented. 



FIG. 2: The renormalized spin-dependent electron energy 
spectra as a function of magnetic field at (a) T=l K and 
(b) r=20 K, with fixed zo^2 nm, d=16 nm, a;=0.002. The 
spin-up (spin-down) levels are denoted by solid (dashed) lines. 



A. Strong lateral confinement 

First we consider a small DMS QD with strong ver- 
tical confinement 2:0=2 nm and small lateral diameter 
d~l6 nm, such that the vertical and in-plane orbital 
energy separations are ~3 eV and ~12 meV, respec- 
tively. The large vertical orbital energy separation en- 
sures only the lowest bound state is relevant to the spin 
relaxation, while the in-plane orbital energy separation 
which is much larger than A^ ensures that the spin-flip 
transitions between the lowest Zeeman doublet usually 
dominate the spin relaxation process. However, if the 
Mn concentration is fairly high and the temperature is 
sufficiently low, the exchange splitting Asd may eventu- 
ally become comparable with the orbital energy separa- 
tion, then the lowest spin-up level may cross spin-down 
excited levels, opening up new spin relaxation channels 
for the spin-up state. 



1. Low Mn concentration 

The Mn concentration is taken as x=0.002, i.e., we 
take the saturated exchange splitting (Asd) sat ~1 meV. 
The renormalized spin-dependent energy spectra for the 
electron at T=l K and T=20 K are shown in Fig. 2. 
At low temperature T=l K (Fig. 2(a)), the thermal- 
averaged Mn spin KS^)! grows rapidly with increasing 
magnetic field (cf. Eq. Q). As a result, the exchange 
splitting Asd increases rapidly to its maximum (~1 meV) 
and saturates, while A^ oc i? increases smoothly, such 
that the electron spin splitting Aq = Asd — Ai reaches 
its maximum at a critical magnetic field, decreases when 
the magnetic field grows stronger, and eventually changes 
its sign at a strong enough magnetic field. In the high 
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FIG. 3: Spin-flip scattering rate of the lowest spin-up level 
and spin-down level as a function of magnetic field at T=l K 
and T=20 K, with the same structure as in Fig. 2. (a) Spin- 
up level, T=l K; (b) spin-down level, T=l K; (c) spin-up 
level, T=20 K; (d) spin-down level, T=20K. The total spin- 
flip scattering rate, the contribution from PZ and DP are 
denoted by solid lines, dashed lines, and short-dashed lines, 
respectively. 



temperature case (see Fig. 2(b)), \{S^\ increases very 
slowly and Aj always dominates, leading to a negative 
Ao whose magnitude increases with increasing magnetic 
field or temperature. 

The SFR's of the lowest Zeeman doublet l/r± are 
shown in Fig. 3 fimction of magnetic field. The 
contributions from piezoelectric coupling (PZ) and defor- 
mation potential interaction (DP) are also indicated by 
the dashed and short-dashed lines, respectively. First, it 
is interesting to notice that 1/t_ in Fig. 3(b) is signifi- 
cantly smaller than l/r-|_ in Fig. 3(a), because the tran- 
sition from (00—) to (OO-h) needs to absorb a phonon, 
but the phonon number is very small at a low tempera- 
ture T—\ K. Second, both l/r-|_ and 1/t_ arc suppressed 
in a strong magnetic field at T=l K, as shown in Fig. 
3(a) and 3(b). The suppression of l/r+ is due to the 
combined effect of small electron spin splitting Aq (cf. 
Fig. 1(d)) and the vanishing phonon absorption factor 
n(|Ao|), while the suppression of 1/t_ is caused by the 
vanishing correlation function G"' (see Fig. 1(a)). We 
also note that the PZ contribution dominates at small 
spin splitting (Ag < 0.3 meV), while the DP contribu- 
tion dominates at large spin splitting (Aq > 0.3 meV), 
which can be clearly seen in Fig. 1(d), Fig. 3(c), 3(d) and 
all subsequent results. This is a direct result of the differ- 
ence between the dependence of the PZ and DP coupling 
constant Q!,y(q) on the wave vector: Q:^^(q) oc 
a^^(q) oc y/g, such that the former (latter) dominates 
at small (large) phonon energy /kjq(= Aq). Finally, the 
electron spin splitting Aq in Fig. 2(a) vanishes at B=0 
T and B «5 T. Correspondingly, the SFR's 1/t± vanish 
in both Fig. 3(a) and Fig. 3(b). 



FIG. 4: Renormalized spin-dependent electron energy spectra 

as a function of magnetic field, at (a) T=\ K and (b) r=50 K, 
with fixed «o=2 nm, d=\% nm, a;=0.01. The solid and dashed 
lines denote the spin-up and spin-down levels, respectively. 



In the high-temperature regime, the spin splitting |Ao| 
increases with increasing magnetic field. Consequently, 
the SFR's 1/t± exhibit the same behaviors, as shown 
in Fig. 3(c) and 3(d). Note, however, the contribution 
from PZ coupling decreases in a strong magnetic field. 
This is caused by the decreasing SV when the electron 
spin splitting Aq exceeds Ad (see Fig. 1(d)). Compared 
with the low-temperature case, wc find that the SFR's 
increase significantly at a higher temperature, due to the 
increasing number of phonons. However, at very strong 
magnetic fields, the absence of high-energy phonon and 
the reduction of , similar to the low-temperature 
case, reduce l/r+ and 1/t_, respectively. 



2. Intermediate Mn concentration 

The Mn concentration is increased to a;=0.01, with 
a saturated exchange splitting {^sd)sat ~5 meV, large 
enough to maintain a positive spin splitting Aq in the 
whole range of the magnetic field i?=0~8 T at low tem- 
perature. Prom the renormalized energy spectra shown 
in Fig. 4, wc sec that As^ (Aj) dominates at T=l K 
(T=50 K) such that Aq is positive (negative) over the 
whole range of the magnetic field. In Fig. 5, the SFR's 
1/t± are plotted as a function of magnetic field for T=l 
K and T=50 K, respectively. At T=l K, when the mag- 
netic field increases, 1/t+ increases to its saturation value 
due to the increase and saturation of Ag^ (and thus Aq), 
while that of the spin-down level is much smaller and 
shows a sharp peak at very weak magnetic field (_B w0.04 
T) and reduces to zero quickly, which is caused primar- 
ily by the small phonon absorption factor n(|Ao|) at low 
temperature and partly by the reduction of the correla- 
tion function G"' in a strong magnetic field. The fur- 
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FIG. 5: Spin-flip scattering rate of the lowest Zeeman doublet 
versus the magnetic field at T=l K and T=50 K, with the 
same structure as in Fig. 4. (a) Spin-up level, T=l K; (b) 
spin-down level, T—1 K; (c) spin-up level, r=50 K; (d) spin- 
down level, T=50 K. The total spin-flip scattering rate, the 
contribution from PZ and DP are denoted by solid, dashed 
and short-dashed lines, respectively. 



FIG. 6: Renormalized spin-dependent electron energy spectra 
vs. the magnetic field at (a) T=l K and (b) r=50 K, with 
fixed zo=2 nm, d—16 nm, a;=0.05. Spin-up and spin-down 
levels are denoted by solid and dashed lines, respectively. 



ther increase of 1/t+ with increasing magnetic field in 
Fig. 5(a) is due to the decrease of the orbital excitation 
energy (in Eq. @, the denominator Ai-^i^ consists of 
two parts, the orbital excitation energy ej^ — ei^ and the 
spin splitting Aq). At T=50 K, 1/t+ and 1/t_ both in- 
crease with increasing magnetic field, due to the increase 
of the spin splitting Aq and the number of phonons. Note 
in Fig. 5(c) and 5(d), the zero-field SFR's do not van- 
ish. This can be understood because the electron can 
transit to excited orbital levels at T=50 K, although the 
spin-fiip transition to the ground orbital level is prohib- 
ited due to vanishing spin splitting Aq. Furthermore, 
we notice the SFR in the case of intermediate Mn con- 
centration (a;=0.01) is several times larger than that for 
low Mn concentration (x=0.002), since 1/t± cx x through 
Eq. {TI) and Eq. ||SJ (note, however, the effect of the Mn 
concentration x is also manifested through changing the 
electron spin splitting Aq). Finally, the decrease of the 
PZ contribution at large electron spin splitting Aq, and 
the resulted crossing of the PZ and DP contributions at 
Aq ~ 0.3 meV, as discussed in the previous subsection, 
is again observed. 



3. High Mn concentration 

In this subsection, the Mn concentration is increased 
further to x=0.05, with a saturated exchange splitting 
{Asd)sat ~16 meV, such that A^d is comparable with 
the in-plane orbital level separation (~ 12 meV). The 
renormalized energy spectra at T=l K and T=50 K are 
shown in Fig. 6. We see from the left panel (T=l K) that 
the spin-up ground orbital level crosses the spin-down ex- 



cited orbital level at a critical magnetic field Be ^2.2 T. 
At a higher temperature T—50 K, Agd is still the dom- 
inant contribution to Aq, but its magnitude decreases, 
such that the energy levels do not cross (see Fig. 6(b)). 

The SFR's of the lowest Zeeman doublet are shown 
in Fig. 7. First we note the SFR's in panel (a), (c), 
(d) are of the same order of inverse nanoseconds, while 
that in panel (b) is much smaller, due to the absence 
of high-energy phonons at low temperature. The most 
significant feature is the sharp peak around the critical 
magnetic field Be in Fig. 7 (a), corresponding to the 
crossing of the spin-up ground level with the first ex- 
cited spin-down level (see Fig. 6(a)). This is because 
the level crossing leads to a resonance in i^oo.oo and, as 
a result, in (see the discussion in the beginning of 

section III). Note, however, 1/t_ in Fig. 7(b) doesn't 
show this resonant behavior, because the resonance of 
-K^oo,oo is suppressed by the vanishing phonon absorption 
factor n(|Ao|). For B < Be, the transition to the lowest 
spin-down level gives the dominant contribution to 1/t_|_, 
which reaches its maximum at B ^0.8 T and decreases 
at stronger magnetic fields (cf. Fig. 1(d)). For B > Be, 
the electron in the lowest spin-up level can transit into 
the first excited spin-down level, opening up a second 
spin-flip channel, and it is just the contribution from this 
channel that dominates in the B > Be regime. This sec- 
ond contribution reaches its maximum value at B w4 T 
and then decreases, which can also be interpreted via Fig. 
1(d). For the T—50 K case, 1/Tf and 1/t_ both increase 
with increasing magnetic field, showing a broad peak at 
B «2.5 T. The peak comes from the competing effect 
of increasing spin splitting Aq (which leads to increasing 
1/t±) against decreasing correlation function G ^,G^ 
and phonon emission (absorption) factor. 
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FIG. 7: Spin-flip scattering rate of the lowest Zccrnan doublet 
as a function of magnetic field at T=l K and r=50 K, with 
the same structure as in Fig. 6. (a) Spin-up level, T=l K; 
(b) spin-down level, T=l K; (c) spin-up level, T=50 K; (d) 
spin-down level, T=50 K. The total spin-flip scattering rate, 
the contribution from PZ and DP are denoted by solid lines, 
dashed lines, and short-dashed lines, respectively. 



B. Weak lateral confinement 

Now wc turn to investigate QD's with weak lateral 
confinement d=40 nm, whose orbital level separation is 
comparable with A^. In this case, with small Mn concen- 
tration or high temperature, A.; makes the main contri- 
bution to the spin splitting. Consequently, the spin split- 
ting Aq is negative, and the lowest spin-down level may 
cross the excited spin-up levels in a strong magnetic field. 
Figure 8 shows the renormalized electron energy spectra 
at T=l K and T—IQ K. In panel (a), where the temper- 
ature is low, the spin-down ground level first crosses the 
spin-up ground level, due to the small exchange splitting 
Asd and low temperature, then it sweeps cross the excited 
spin- up levels, due to the large Aj compared with the or- 
bital excitation energy in a strong magnetic field. When 
the temperature increases to T=10 K (see Fig. 8(b)), 
the exchange splitting is suppressed and Aj always dom- 
inates. The lowest spin-down level crosses the excited 
spin-up levels but the crossing between the Zeeman split 
doublet doesn't occur. 

Figure 9 shows 1/t± as a function of magnetic field at 
T=l K and T=10 K, respectively. In panel (a) and (b), 
the low-field behaviors of l/r± resemble those of strong 
lateral confinement and low Mn concentration (see Fig. 
3(a) and 3(b)). At higher temperature, 1/t± exhibit 
many peaks at higher magnetic fields (indicated by the 
arrows), which are caused by the aforementioned level 
crossings. However, in very strong magnetic fields, the 
peaks are suppressed by the phonon absorption factor 
(for 1/t+) and the correlation function (for 1/t_). 

At a higher temperature T=10 K, the resonances of the 



FIG. 8: Renormalized spin-dependent electron energy spectra 
vs. the magnetic fleld at (a) T=l K and (b) T=10 K, with 
fixed zo=2 nm, d=40 nm, a;=0.001. Spin-up and spin-down 
levels are denoted by solid and dashed lines, respectively. 




B(T) B(T) 



FIG. 9: Spin-fiip scattering rate of the lowest Zeeman doublet 
as a function of magnetic field at T=l K and T=10 K, with 
the same structure as in Fig. 8. (a) Spin-up level, T=l K; 
(b) spin-down level, T=l K; (c) spin-up level, T=10 K; (d) 
spin-down level, T=10 K. The total spin-flip scattering rate, 
the contribution from PZ and DP are denoted by solid lines, 
dashed lines, and short-dashed lines, respectively. 



kernel -K^oo.oo are less siippresscd and more pronounced 
peaks arises in 1/t±, leading to short spin lifetimes of 
the order of nanoseconds, compared with the microsec- 
ond scale in the T=l K case. 



C. Temperature effect 

In the above, we have observed that the temperature 
plays an important role in determining the SFR through 
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Temperature T (K) Temperature T (K) 



FIG. 10: (a) Renormalized energy of the lowest spin-up (solid 
line) and spin-down (dashed line) levels vs. the temperature, 
(b) and (c) show the the spin-flip scattering rate of the spin- up 
and spin-down levels vs. the temperature, respectively. The 
total spin-flip scattering rate, the PZ and DP contribution are 
denoted by solid, dashed and short-dashed lines, respectively. 
Here we take 3—4 T, zo=2 nm, d—16 nm and x=0.01 in the 
numerical calculation. 



changing the electron spin splitting Aq, the correlation 
function G"' , G ^ and the phonon emission (absorp- 
tion) factor. Taking a small QD (2:0=2 nm, d—16 nm) 
for example, we plot the renormalized energy spectrum 
and 1/t± as a function of temperature at B=A T in Fig. 
10. It is interesting to notice that the Zeeman split dou- 
blet crosses each other at an elevated temperature (see 
Fig. 10(a)), due to the reduction of the exchange split- 
ting Asd- From Fig. 1(d), we see that the quantity ST 
decreases with decreasing Ao for small Aq. This effect, 
together with the phonon emission factor, which shows a 
sharp peak at Ao=0, results in the non-monotonic tem- 
perature dependence of l/r+ shown in Fig. 10(b). In Fig. 
10(c), the low-temperature SFR for the spin-down level 
1 /t_ vanishes due to the absence of high-energy phonons 
and the vanishing correlation function . Note at 
T «34 K, the SFR of both spin-up and spin-down levels 
vanishes, due to the vanishing spin splitting Aq. 



D. Dependence of the SFR on the confinement 

Both the vertical and lateral confinement of the QD 
can affect the spin relaxation significantly through vary- 
ing the electron wave function. The effect of vertical 
confinement on 1/t± comes from the form factor Z{q) 
(see Eq. (IA.5|l ) and the overlap integral (see Eq. (|A.2|l ). 
It can be seen from Fig. 1(c) that the quantity ST is 
roughly proportional to 1/2:0, such that the SFR's should 
also show the same behavior, which can be seen in Fig. 
11(a) and 11(b). Note here the largest zq (20nm) still 



FIG. 11: (a) and (b) show l/r± vs. the vertical conflning 
length 20 at fixed B=l T, T^l K, d=16 nm, a;=0.001. (c) 
and (d) show l/r± vs. the lateral diameter d at fixed _B=4 
T, r=l K, x=0.002. The total spin-fiip scattering rate, the 
PZ and DP contributions are denoted by solid, dashed and 
short-dashed lines, respectively. 



sustains a vertical orbital level separation of ^ 30 meV, 
such that the influence of higher subbands on the spin 
relaxation is negligibly small. The approximate relation- 
ship 1/t± oc 1/20 comes from two factors. First, the 
s-d exchange scattering amplitude with one Mn ion is 
proportional to I/zq which, when squared (in the Fermi 
golden rule) and averaged over all the Mn sites, leads to 
the 1/20 dependence of the spin-flip scattering rate to a 
given final state. Second, the spin- flip channel (i.e., the 
number of final states) doesn't increase provided zq is 
small enough such that only the lowest bound state is 
relevant. We notice that G. Bastard et al. performed a 
theoretical calculation of the SFR of subbands in DMS 
quantum wells, and similar dependence of the SFR on 
the well width is predicted. 

The effect of the lateral confinement strength, charac- 
terized by the lateral diameter d, on the spin relaxation 
is shown in Fig. 11(c) and 11(d). The dependence 
of 1/t± can be appreciated as follows. From the Ap- 
pendix and the definition d = 2-^/ h/ {m*ujQ), we see the 
orbital level separation S is roughly proportional to 1/d^, 
while 1/t± is inversely proportional to (5^ (see Eq. Q) 
when the spin splitting is small, so we expect that 1/t± 
should be approximately proportional to d"*, although the 
precise dependence of l/r± on c? is also affected by the 
phonon-induced transition rate T (cf. Eq. (|A.3p . (|A.4p ). 

The dependence of the SFR on the QD size l/r± cx 
ZQ^d'^ caused by phonon-mediated s-d exchange scatter- 
ing is quite different from those caused by other spin 
relaxation mechanisms in nonmagnetic semiconductor 
Qj),gi^,2ij22 -y^g note that this relationship can be de- 
duced from the work by Nazarovf^ where the spin relax- 
ation in nonmagnetic semiconductor QD's is considered, 
but the magnitude of the SFR in our results is several 
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orders of magnitude higher. 



jor State Basic Research Project No. G001CB3095 of 
China. 



IV. CONCLUSIONS 



Based on second-order perturbation theory, we have 
investigated the SFR caused by the phonon-mediated 5- 
d exchange interaction of the lowest Zeeman spht doublet 
in II- VI DMS QD's as a function of magnetic field, as well 
as the dependence of the SFR on the Mn concentration, 
dot size and temperature. We found the SFR increases 
with increasing Mn concentration and electron spin split- 
ting Aq. Increasing the lateral dot size leads to enhanced 
SFR while increasing the vertical dot size reduces the 
SFR for a small QD. The dependence of the SFR on the 
magnetic field and temperature shows non-monotonic be- 
haviors, due to the competing effect between the electron 
spin splitting, the phonon emission (absorption) factor 
and the correlation function of the Mn ions. It is inter- 
esting to notice that the spin relaxation of both spin- 
up and spin-down electrons is suppressed in the case of 
strong magnetic field and low Mn concentration at low 
temperature. 
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APPENDIX: EXPRESSIONS OF S AND F 



The orbital part of Hq gives the Fock-Darwin states 



KmiP, = 4=e""^ • Rnmip) («, |m| = 0, 1, 2, • • • ), 

V -^TT 

(A.1) 



with corresponding orbital energy Snm = (2n + |m| + 

l)huj + ■mhLUc/2, where Iq = y/h/Jm*Lu), L'^i^-) is the 
generalized Laguerre polynomial. For convenience, we 
also introduce n+ = n + (|m| + m)/2, ri- = n + (|m| — 
m)/2. 

The dimensionless overlap integral is 



_ . / ni!n2!n3!n4! fl^ 

/ e-2-(V5)'™^i+i'"=i+i'"^i+i'"*i47ii(x)47i(x)47i(x)Li7; 

^0 



(A.2) 



where ^ = zq J dz \x{^)t — 3/2. The phonon transition is 
rate due to piezoelectric coupling to the acoustic phonon 



J 



I (ni-K,<)!(ni-,<)!(n2+,<)!(n2-,<)! , J„,^_„' | + |„,__„M-Ff (fM^ /a Q^ 

i-m2,m[-m'\ 1 TTT TTT TTT ^1^^^) ' ' ' V^-^^ 

1 ^ >' (ni+,>)!(ni_,>)!(n2+,>)!(n2-,>)' 



Airhp 



JV-l-l 



de. 



The contribution from the deformation potential interac- 



tion is 
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I1L1I2I2 



'(ni+,<)!(ni_,<)!(n2+,<)!(n2_,<)! 



^mi—m2,mi- 



(ni+,>)!(ni_ 



,>)!("2+,>)!('^2-,> 
,272 „;_2/) ^272 



(-1)1 



"2+- 



2+| + h2- 



(A.4) 



/ (sin^)^+i |Z(giCos^)| exp(-^^-^- 

^0 z 4 



In the above, we have used the short notation I to denote 
the quantum number {n,m) and nj,< (nj,>) to denote 
min {nj,n^} (max jrij, n^}), e.g., stands for (ni,TOi) 
and ni+,> stands for max {ni+, n'j^^}. Other quanti- 



ties are N = |ni_|_ — 
n2- — n 



\ni- 



T-i-l + P2+ - "-2+I + 
= \ Aii> I /{hC,y), C„ {v = is the longi- 
tudinal or transverse sound velocity, 6 is the polar angle 
of the phonon vector q, A,, [0) is the anisotropy function 
of the piezoelectric interaction, with the dependence on 
the the azimuth angle <^ averaged out. 



I' 



,iqz 



AirHie'i^" - 1) 
qzo [{qzor - (27r)2 



(A.5) 



the 

|ni+-"i+| 



form 



and S(i is the deformation potential constant. 



factor, 

' (.'r)Z/i2+,< 



"2- 



For the spin-flip transitions between the lowest Zeeman 
doublet |00-f) and |00— ) , the following overlap integral 
and phonon transition rates are used: 



c _j (m +n2 |mi|)! 

'J00,nimi,00,n2TO2 — ''mi,m2 „ , ,, , , ; 77 \ nrrT i rrr n ;2 ' V-^-'-'J 

2"i+"2+l'»ilv'ni!n2!(ni -h |mi|)!(n2 -h |m2|)! 27rzo^g 

pPX _ s: (_l)ni-Fn2 (e/i.i4)^ Jog;. . 2(ni-Hn2 + |mi |) g'^ / a 7^ 

- V.,^2 ^ |^,|)!n2!(n2 + |m2|)! 47r?lp 2 ^ ^""-'^ 

rA,(g)|Z(g.cos^)|^exp(- ^°^'!^"'^ )^( ^°'^'f"'^ )(sin^)^("^+"-+l-"^l)+^dg, 
Jo 2 4 



r-DP _X ^) ^d^l c g'^0 \2(ni+n2-Hmi|) /aox 

J- 00,nimi,00,n2m2 — ''TOi,m2 y j 77 — j p-77 — j m zlTr/i V 9 ^ V-^-°/' 

A/ni!n2!(ni + \mi\)\(n2 + |m2|)! 47rftpC; / 



(sin l)+i I ^(5; COS 0) I ^ exp(- ^i^i^)^( ^ 



and >S'nimi,ra2m2,oo,oO) roofniTni,n2m2,oo; roc^nimi,n2m2,oo -'^oo^imi,oo,n2m2 ' respectively, by replacing ^^j^^^^ with 

can be obtained from 5oo,nimi,00,n2m2> ^Eo,nirmfiO,n.2m2-: ^^i.-ms- 
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